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ABSTRACT 

We compare evolutionary models for protoplanetary discs that include disc winds with 
observational determinations of the disc lifetime and accretion rate in Taurus. Using 
updated estimates for stellar ages in Taurus, together with published classifications, we 
show that the evolution of the disc fraction with stellar age is similar to that derived 
for ensembles of stars within young clusters. Around 30 percent of stars lose their 
discs within 1 Myr, while the remainder have disc lifetimes that are typically in the 
1-10 Myr range. We show that the latter range of ages is consistent with theoretical 
models for disc evolution, provided that there is a dispersion of around 0.5 in the log 
of the initial disc mass. The same range of initial conditions brackets the observed 
variation in the accretion rate of Classical T Tauri stars at a given age. We discuss the 
expected lifetime of discs in close binary systems, and show that our models predict 
that the disc lifetime is almost constant for separations d £ 10 au. This implies a low 
predicted fraction of binaries that pair a Classical T Tauri star with a Weak-lined T 
Tauri star, and is in better agreement with observations of the disc lifetime in binaries 
than disc models that do not include disc mass loss in a wind. 

Key words: accretion, accretion discs — stars: pre-main-sequence — planetary 
systems: protoplanetary discs — open clusters and associations: general 



1 INTRODUCTION 

Observations of young stars in star forming regions such 
as Taurus Aurigae show that there is a wide dispersion in 
the lifetimes of circumstellar discs, in the range 1-10 Myr 
(Strom et al. 1989; Haisch, Lada & Lada 2001). Although 
disc lifetimes of this order are consistent with the meteoretic 
evidence from our own Solar System (Russell et al. 1996) , the 
reason for the large dispersion in the disc lifetime remains 
unknown. The obvious controlling parameter would be the 
stellar mass, but in fact observations of low mass stars in 
young clusters show little correlation between mass and disc 
lifetime (e.g. Rebull et al. 2002). 

The simplest disc models also fail to reproduce other 
properties of circumstellar discs. One problem is to match 
the rapid dispersal of the disc. An accretion disc evolves 
on a single, viscous timescale, which in order to reproduce 
the aforementioned disc lifetimes needs to be of the order of 
10 6 yr. Dispersal of the disc on this timescale would lead to 
the prediction of a substantial class of objects with proper- 
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ties intermediate between accreting Classical T Tauri stars 
(CTTS) and discless Weak-lined T Tauri stars (WTTS) (Ar- 
mitage, Clarke & Tout 1999). In fact, few stars are caught 
in the act of clearing their discs, which constrains the disc 
dispersal time to be short - of the order of 10 5 yr (Simon & 
Prato 1995; Wolk & Walter 1996; Duvert et al. 2000). 

A second problem is that disc evolution is largely in- 
different to the presence or absence of close binary compan- 
ions. The viscous timescale increases with the characteris- 
tic radius of the disc, so discs hemmed in by close binary 
companions ought to evolve more rapidly than discs around 
relatively isolated stars. Although examples are known of bi- 
naries that appear to have dissipated their discs at an early 
epoch (Meyer et al. 1997) , surveys suggest that binary com- 
panions as close as 10 au generally have little effect on the 
lifetime or mass accretion rate of the disc surrounding the 
primary (Simon & Prato 1995; Bouvier, Rigaut & Nadeau 
1997; Chez, White & Simon 1997; White & Chez 2001). 

Recently, Clarke, Gendrin & Sotomayor (2001) have 
shown that models which include mass loss from the disc 
lead to a dispersal timescale that is much shorter than the 
viscous evolution timescale, reproducing the two-timescale 
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behavior required by the observations. In this paper we 
demonstrate that the same class of models, when combined 
with a modest range of initial disc parameters, are quantita- 
tively consistent with the wide dispersion in disc lifetimes. 
In Section 2 we set out the disc model and show that when 
mass loss is included, the evolution of the accretion rate de- 
pends mainly on the initial disc mass and viscous timescale 
at large radius, and only weakly on the radial dependence 
of disc parameters such as the surface density and viscos- 
ity. In Section 3, we experiment with the range of initial 
conditions that are required to reproduce the observed age 
distributions of CTTS and WTTS in Taurus- Aurigae, and 
show that this range also reproduces the observed scatter in 
accretion rates in CTTS inferred by Gullbring et. al. (1998). 
In Section 4, we apply the model to binary star systems, and 
show that due to the wind mass loss from large radii in the 
disc, the disc lifetime may increase with decreasing disc bi- 
nary separation over the binary separation range 10 — 100 au. 
Section 5 summarizes our conclusions. 



2 DISC MODEL 

2.1 Evolution equation 

The evolution of the surface density, E(r, t), of a geomet- 
rically thin accretion disc is described by (Lynden-Bell & 
Pringle 1974; Pringle 1981), 
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where the term E w i n d represents the rate of mass loss per 
unit surface area of the disc in the equatorial plane. We 
assume that any mass lost leaves the disc with Keplerian 
specific angular momentum, and thus does not drive radial 
inflow or outflow in the remaining disc gas. For the kinematic 
viscosity v we adopt a simple power-law in radius, 

,/3 



v = Vor 



(2) 



with vq and /3 constants. Within the context of a disc mod- 
els (Shakura & Sunyaev 1973), the viscosity in a protoplane- 
tary disc is determined by the intrinsic efficiency of angular 
momentum transport (i.e. the value of a); by the relative 
contribution of viscous heating and stellar irradiation; and 
by the opacity which controls the vertical disc structure. A 
passive irradiated disc, in which both the central temper- 
ature and the surface temperature decline with radius as 
T c oc r -1 / 2 , has a viscosity (assuming a constant a) which 
scales as v on r. More elaborate calculations of protoplane- 
tary disc structure - though still within the a disc formalism 
- likewise yield /3 ~ 1 (e.g. Bell et al. 1997; D'Alessio et al. 
1998). Only limited observational constraints are available at 
the (small) radii of interest, though j3 = 1 is consistent with 
VLA observations of TW Hydrae by Wilner et al. (2000). 

In a steady-state, the surface density profile correspond- 
ing to a given mass accretion rate M is, 



M 
3^ 
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where we have assumed standard zero-torque boundary con- 
ditions at n n , the inner edge of the disc. 

Equation Q is a diffusion equation for the surface den- 
sity. Provided that the surface density evolution is controlled 



by angular momentum transport rather than mass loss, the 
characteristic timescale for evolution is the viscous time, 

r 2 

t„=-. (4) 

Since the viscous timescale increases with radius, the evo- 
lution of a disc described by equation Q speeds up at late 
times if the outer edge of the disc is truncated by mass loss. 

Standard numerical methods (e.g. Pringle, Verbunt & 
Wade 1986) are employed to solve equation Q on a non- 
uniform radial mesh. We use zero-torque (E(rj n ) = 0) 
boundary conditions at the inner edge of the disc. The spa- 
tial domain is large enough that the outer boundary condi- 
tion has no influence upon the results. Formally, we set the 
radial velocity v r = at this boundary. 



2.2 Disc mass loss 

Observations provide direct evidence for the existence of two 
types of disc outflow, jets and photoevaporative flows. Of 
these, photoevaporative flows, which lead to mass loss from 
the outer regions of the disc, are the most obviously promis- 
ing for producing two-timescale disc evolution and rapid 
disc dispersal (Clarke, Gendrin & Sotomayor 2001; Mat- 
suyama, Johnstone & Hartmann 2002). We consider mass 
loss via photoevaporation exclusively in this paper, while 
noting that almost identical behaviour is also possible in 
models of magnetically driven outflows in which mass loss 
extends beyond the inner few au of the disc (e.g. Ostriker 
1997). 

Photoevaporative flows are thermal winds that occur 
when the disc surface is heated by exposure to ultraviolet 
radiation (Hollenbach et al. 1994; Johnstone, Hollenbach & 
Bally 1998; Hollenbach, Yorke & Johnstone 2000; Richling 
& Yorke 2000). The heated gas can escape at radii r > r g , 
where the critical radius r g is given (approximately) by the 
condition that the sound speed c 3 of the hot gas exceeds the 
local escape velocity, 

r g ~ — g-. (5) 

For parameters appropriate to T Tauri stars, r g ~ 5— 10 au. 
The flux of photoevaporating radiation can be external, as 
in Orion where UV radiation from massive stars is responsi- 
ble (Churchwell et al. 1997; Johnstone, Hollenbach & Bally 
1998; Scally & Clarke 2001; Henney et al. 2002). Alterna- 
tively, a lower level of mass loss can result from the UV 
and X-ray flux of the low mass star itself (Shu, Johnstone 
& Hollenbach 1993). 

In Taurus, there are no massive stars which could gen- 
erate a strong external flux of photoevaporating radiation. 
We therefore adopt a simple mass loss prescription appro- 
priate for the case where the mass loss is driven by local 
irradiation (Hollenbach et al. 1994). The mass loss per unit 
area of the disc, £ W md, is a function of radius but not of 
time , and is given by, 



0, 



r <r„ 



1 Implicitly, this assumes that the photoevaporating flux is not 
tied to the accretion rate. Models by Matsuyama, Johnstone & 
Hartmann (2002) show that photoevaporation driven solely by the 
star-disc accretion shock fails to disperse the disc quickly enough. 
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r / au 

Figure 1. The accretion rate (upper panel) and the surface den- 
sity profile (lower panel, plotted at t = 0.5 Myr) for disc models 
with /3 = 3/2 (solid line), (3 = 1 (short dashes), and ,3 = 1/2 
(long dashes). The other parameters of the models are as given 
in Table 1. Models with very different viscosity laws and surface 
density profiles can yield almost identical accretion rate histories. 

S„md oc r~ 5/2 , r > r g . (6) 

We express the normalization of the mass loss via the pa- 
rameter M w i n d, which is defined as the total mass loss rate if 
the disc extends to 25 au. The instantaneous mass loss rate 
will therefore differ from this depending upon the outer disc 
radius. We take r g — 5 au. 

2.3 Model parameters 

The model described above is denned by two parameters de- 
scribing the disc physics (yo and 0), two describing the wind 
(M w i n d and r g ), and two for the initial conditions (iWjnit and 
the initial disc mass) 2 . Sensible values for most of these pa- 
rameters can be estimated from existing observations. We 
adopt Afmit = 5 x 10 -8 Meyr -1 , since most T Tauri stars 
have accretion rates of this order or lower (Gullbring et al. 
1998). With this initial accretion rate, a disc lifetime of 1- 
10 Myr mandates an initial disc mass of « O.IMq. After 
2 Myr of evolution, this initial disc mass has declined (see 
Table 1) to a few hundredths of a Solar mass, a value which 
is broadly consistent with estimates of T Tauri disc masses 
(Beckwith & Sargent 1991). The assumption that mass loss 
is driven by photoevaporation fixes r g to be 5-10 au, and 
the requirement that this mass loss drives rapid dispersal of 
the disc suggests that M w i n d be no more than a couple of 

2 Equivalently, we could parameterize the initial conditions in 
terms of the initial mass and angular momentum content of the 
disc. 



orders of magnitude smaller than Mintt (Clarke, Gendrin & 
Sotomayor 2001). 

Determining the parameters for the disc viscosity is 
more difficult. We have already noted that theoretical mod- 
els for protoplanetary discs favour values for /3 around unity, 
but no-one would claim that this parameter is as yet accu- 
rately determined. We have therefore chosen for our fiducial 
models combinations of Vo and /3, shown in Table 1 along 
with the other parameters of the disc model, that yield a 
disc lifetime of approximately 5 Myr. This choice of disc 
lifetime fixes the viscous time at a radius of around 10 au, 
but does not constrain either vo or /3 individually. Indeed, 
as shown in Figure 1, disc models that include a disc wind 
have only a limited period during which the accretion rate 
drops as a power-law denned by the value of (3 in the vis- 
cosity law. As a consequence of this, models with different 
choices of (3 can be constructed that yield almost identical 
evolution of the disc accretion rate, even though they have 
entirely different surface density profiles. This degeneracy 
implies that it is not possible to discriminate observation- 
ally - even if we had perfect data - between models with 
different /3, using only measurements of the evolution of the 
accretion rate with time. We cannot, therefore, use the ac- 
cretion rate evolution measured by Gullbring et al. (1998) to 
constrain our models in the same way as was done for disc 
models without mass loss by Hartmann et al. (1998) and 
Stepinski (1998). In the following discussion, we will leave j3 
as an unconstrained parameter, and consider the three mod- 
els outlined in Table 1 as equally valid possible descriptions 
of the disc. 

To illustrate how the results depend upon some of the 
other model parameters, we show in Fig. [5] how the accre- 
tion rate evolves with time for different choices of the initial 
disc mass, initial accretion rate, and disc mass loss rate. We 
start with a disc of mass Mdiac, and accretion rate Minit- 
The disc has an initial surface density profile defined by 
equation out to the radius where the enclosed mass is 
Mdisc, and E = at greater radii. As is obvious from Fig. [21 
the two-timescale behaviour discussed by Clarke, Gendrin & 
Sotomayor (2001) occurs in a range of disc models which in- 
corporate mass loss. The decline of the accretion rate steep- 
ens dramatically in all models at the epoch when mass loss 
in the disc wind becomes important. 



3 THE DISC FRACTION IN TAURUS 

In a recent study, Haisch, Lada & Lada (2001) determined 
(using JHKL band photometric surveys)the fraction of stars 
possessing discs as a function of the previously published 
mean cluster age. They found an evolution in the disc frac- 
tion shown in the lower panel of Figure 4. However, since 
the spread in stellar ages can be comparable to the mean 
age for clusters that are only ~ 10 6 yr old, it is not imme- 
diately obvious whether the smooth decline in disc fraction 
with 'age' seen in the Haisch, Lada & Lada (2001) sample 
reflects primarily, 

• A true dispersion in the lifetime of circumstellar discs. 

• A convolution of a constant disc lifetime with an ex- 
tended epoch of star formation in each cluster. 
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/3 = l/2 


/9 = 1 


,8 = 3/2 


v(r = 1 au) 


6.5 x 10 13 cm 2 s _1 


7.5 x 10 13 cm 2 s- 1 


1.3 X 10 14 cm 2 s" 1 


Initial disc mass 


0.15 Mq 


0.125 M 


0.10 M 


Disc mass at 2 Myr 


0.07 Mq 


0.045 M Q 


0.027 Mq 


Viscous time at 10 au 


3.9 Myr 


4.8 Myr 


3.9 Myr 


Inner radius for disc mass loss 


5 au 


5 au 


5 au 


Afinit 


5 X 10" 8 Mgyr" 1 


5 X 10~ 8 Mgyr" 1 


5 X 10" 8 Mgyr" 1 


-^wind 


5 X 10" 9 Mgyr" 1 


5 X 10~ 9 Mgyr" 1 


5 X 10~ 9 Mgyr" 1 



Table 1. Summary of parameters for the three models shown in Fig. 1. All of these models are broadly consistent with the observed 
evolution of the accretion rate and disc mass in young clusters such as Taurus. 




t / yr t / yr t / yr 



Figure 2. Illustration of how the accretion rate evolution depends upon the model parameters. The fiducial model has v <x r 3 / 2 , 
an initial disc mass M disc = 0.1 Mq, M init = 5 X 10~ 8 Mgyr^ 1 and M wind = 5 X 10~ 9 Mgyr^ 1 . We adopt r in = 0.066 au, and 
r out = 100 au. The leftmost panel shows models with mass loss rates M wind = 5 X 10 _11 ,5 X 10 _10 ,5 X 10 _9 ,5 X 10 -8 Mgyr -1 . 
The centre panel shows the effect of changing the initial disc mass to 0.025,0.05,0.1,0.2 Mq. The rightmost panel shows models with 
Mnit = 10~ 8 ,2.5 X 10" 8 ,5 X 10- 8 ,10" 7 Mgyr" 1 . Rapid dispersal of the gas is a generic feature of all of these models. 



To distinguish between these possibilities, it is necessary to 
date the stars within young clusters on a star by star basis. 
We present such an analysis below, and show that in Tau- 
rus the evidence suggests that for low mass stars there is a 
significant intrinsic dispersion in disc lifetime. 



3.1 Data 

To study the evolution of the disc fraction with stellar age, 
we make use of data compiled for a study of star forma- 
tion in Taurus by Palla & Stahler (2002). The sample of 
151 stars includes approximately 10 stars from Wichmann 
et al. (2000) whose pre-main-sequence status is confirmed 
by lithium observations, and which fall within the bound- 
aries of the 12 CO map that contains the bulk of the Taurus- 
Auriga population. ROSAT All Sky Survey sources that fall 
at larger distances have not been included. Of this sample, 
67 stars are classified as CTTS, and 57 are WTTS. For 27 
stars discovered in more recent surveys we were unable to 



find published classifications, and these are excluded from 
further analysis. 

For each of the stars in the CTTS and WTTS sam- 
ples, we have estimated ages using the evolutionary tracks 
of Palla & Stahler (1999). In general, age estimates for pre- 
main-sequence stars are subject to considerable uncertainty, 
due to different assumptions as to the initial conditions and 
accretion history (Tout, Livio & Bonnell 1999; Hartmann 
2001; Baraffe et al. 2002). In the specific case of Taurus, 
Hartmann (2003) has argued that the ages of the higher 
mass stars (with masses > 1 Mq) are especially uncertain, 
due to different treatments of the stellar birthline (e.g. Hart- 
mann, Cassen & Kenyon 1997; Palla & Stahler 1999) and 
possible foreground contamination. We therefore further re- 
strict the sample under consideration to those stars with 
effective temperatures T e < 4350 K, leaving a final sample 
of 57 CTTS and 41 WTTS. With the effective temperature 
cut, these are all low mass stars. 

Figure[3]shows histograms of the age distribution of the 
low mass CTTS and WTTS in Taurus. The most immedi- 
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Figure 3. Distribution of the ages of stars classified as Classical 
T Tauri stars (upper panel) and Weak-Lined T Tauri stars (lower 
panel) in Taurus. Only stars with effective temperature T e < 
4350 K are included in the histograms. 



ately obvious conclusion to be drawn from the Figure is that 
there is substantial overlap in the derived ages of the CTTS 
and the WTTS. However, there is also evidence that the 
WTTS are significantly older as a class. A KS test reveals 
that the age distributions of the WTTS and CTTS samples 
differ at more than the 2a level (the probability that the 
distributions are the same is approximately 0.03). 

Although the empirical conclusion that the WTTS arc 
significantly older than the CTTS falls short of statistical 
proof, the evidence is consistent with our expectations if 
CTTS turn into WTTS. Indeed, if no such difference were 
found, one would either have to posit no evolutionary link 
between the two (thus begging the question of what CTTS 
should turn into) or else appeal to contrived models in which 
the lifetime of discs were shorter in the case of stars born 
more recently. We shall show below that the Taurus data is 
instead compatible with the simple hypothesis that CTTS 
turn into WTTS, but that the clock for this process varies 
from star to star. In Chamaeleon I, by contrast, it has been 
claimed that the age distribution of CTTS and WTTS is 
indistinguishable (Lawson, Feigelson & Huenemorder 1996). 
It is now known, however, that this dataset was missing 
a number of young CTTS (Persi et al 2000), and a more 
complete census of this region may not ultimately yield an 
answer different from the Taurus result reported here. 

Figure [I] shows the CTTS fraction in Taurus as a func- 
tion of stellar age. We have binned the data into 5 bins, 
evenly spaced in the log of the estimated stellar age. Most 
of the transition between CTTS and WTTS occurs for ages 
between 1 Myr and 10 Myr. At the youngest inferred ages, 
CTTS are dominant, but there is still a significant admix- 
ture of WTTS. We note, however, that for such young ages 
- less than a Myr - theoretical uncertainties in the age are 



Figure 4. The disc fraction as a function of age in Taurus and 
other young clusters. The upper panel shows the Classical T Tauri 
fraction /ctts = -Wctts /C^CTTS + -^wtts) as a function of 
stellar age, using published classifications and our updated es- 
timates for the ages. The lower panel compares the Taurus data 
(triangles, shifted slightly to the right for clarity) to the evolution 
of the disc fraction measured by Haisch, Lada & Lada (2001) for 
young clusters (points). 



probably comparable to the age itself (Baraffe et al. 2002). 
Any conclusions we might draw as to the evolution of the 
CTTS fraction prior to 1 Myr should therefore be treated 
with caution. 

Figure 0] also shows the evolution of the disc fraction as 
measured by Haisch, Lada & Lada (2001) in a sample of 6 
young clusters with estimated ages ranging between 0.3 Myr 
and 30 Myr. Within the errors, which are substantial (arising 
both from the limited number of stars, and from the uncer- 
tainty in cluster age which is not shown in the Figure), the 
same trend with age is seen as for the Taurus data. 



3.2 A model for the dispersion in disc lifetime 

A dispersion in the disc lifetime could be caused by several 
factors, including different disc initial conditions or varying 
rates of mass loss. For photoevaporation, the mass loss rate 
is predicted to vary only weakly (as the square root) with the 
ionizing flux. We have therefore explored whether the data 
shown in Fig.^Jcan instead be explained as arising from dif- 
ferent initial conditions for the discs. We focus on the effect 
of different initial disc masses, as Fig. [5] indicates that the 
disc lifetime changes markedly as we vary this parameter. 

Fig. |S] shows how the disc lifetime varies with the ini- 
tial disc mass. We define the disc lifetime as the time at 
which the accretion rate first drops below 10~ 10 Moyr -1 
(once this epoch is reached, as shown in Figure 1, mass loss 
in the wind dominates the evolution, and the accretion rate 
is dropping precipitously). Apart from changing Maine, the 
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Figure 5. The disc lifetime as a function of the initial disc mass, 
for models with (3 = 3/2 (solid line), (3 = 1 (short dashes), and 
(3 = 1/2 (long dashes). The initial accretion rate is assumed con- 
stant for all models. 
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Figure 6. Evolution of the disc fraction predicted by models with 
(3 = 3/2 (upper solid line), (3 = 1 (short dashes), and (3=1/2 
(long dashes), assuming in each case that there is a dispersion of 
0.5 in the log of the initial disc mass. Better, but ad hoc fits to 
the data (lower curves), are obtained by assuming that « 30% of 
stars lose their discs via other processes within the first Myr. 



other parameters are kept fixed at the values shown in Ta- 
ble 1 (i.e. we change the initial conditions for the disc, but 
keep the parameters describing the disc viscosity and disc 
wind constant). For all three disc models, the disc lifetime 
has an approximately power-law dependence upon the initial 
disc mass, with an index of w 0.5. 

Figure |S| shows the predicted disc fraction as a function 
of age, under the assumption that log(A/disc) is distributed 
as a Gaussian around the central values given in Table 1. 
We find that the general trend in the plot of disc fraction 
with age can be reproduced provided that there is a 1 sigma 
dispersion of approximately 0.5 in the log of the initial disc 
mass. A dispersion in initial disc mass of this magnitude 
leads to almost all stars losing their discs between 1 Myr 
and 10 Myr, as required to fit the evolution of the CTTS 
fraction in Taurus, and the evolution of the disc fraction in 
the clusters studied by Haisch, Lada & Lada (2001). 

Although the models do a reasonable job in reproduc- 
ing the width of the transition seen in Fig. |SJ they fail in 
that they predict a 100% disc fraction for the youngest 
stars, those with ages less than a Myr. Indeed, from Fig. 
we find that a disc lifetime as short as (say) 0.5 Myr re- 
quires an initial disc mass that is very small - of the order 
of few xlO -3 Mq. As noted previously, one possibility is 
simply that the ages of the youngest WTTS have been un- 
derestimated. We find, however, that the spatial distribution 
of these young WTTS, relative to the cluster gas, is more 
consistent with the distribution of the young CTTS (Palla 
& Stahler 2002) than of the older stars. Another possibil- 
ity is that these objects, classified as WTTS on the grounds 
of their equivalent width in H a , have highly variable line 
emission, and have thus been classed as WTTS even though 
they are not bona fide discless objects. However, we find 
that these systems are also lacking significant excesses both 
in the near and far infrared, thus confirming their discless 
status. These considerations support the contention that the 
apparently young, discless stars really are young and disc- 
less, and constitute a separate population which have lost 
their discs at a earlier epoch due to processes that do not 
affect the majority of stars. An ad hoc model in which 30% 
of stars lose their discs within a Myr, while the rest evolve 
in the way which we have modelled, is shown in Fig. |S| and 
is consistent with the available observations. 

If stars form in isolation, then it is hard to understand 
why a significant fraction of them should lose their discs 
very early. Rapid early accretion onto the star is expected 
to cease when the disc first becomes gravitationally stable, 
at a mass of the order of 0.1 M,, after which time slower 'vis- 
cous' evolution continues. Within this framework, initial disc 
masses that are smaller than average by one to two orders 
of magnitude are puzzling. In an environment where stars 
are strongly sub-clustered at birth (Scally & Clarke 2002), 
however, close encounters between stars can efficiently re- 
duce the disc mass, both directly by unbinding disc mate- 
rial (Clarke & Pringle 1993), and indirectly by reducing the 
outer disc radius and hastening viscous evolution (Armitage 
& Clarke 1997). Recent simulations of clustered star forma- 
tion (Bate, Bonnell & Bromm 2003) show that the required 
interactions may well be extremely common. In the Bate et 
al. (2003) calculation, most stars have encounters with an- 
other star with a minimum separation distance of a few au 
or smaller. The young WTTS seen in Taurus may simply be 
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Figure 7. The accretion rate of Classical T Tauri stars in Taurus, 
as estimated by Hartmann et al. (1998), plotted against the stellar 
age. The heavy central curve shows the predicted evolution of the 
accretion rate for the ft = 3/2 disc model outlined in Table 1. 
The upper and lower curves show the accretion rate for models 
in which the initial disc mass has been increased or decreased by 
a factor of three. 



Figure 8. The predicted disc lifetime in binary systems, for mod- 
els with (3 = 3/2 (filled triangles), = 1 (open squares), and 
(3 = 1/2 (crosses). We assume that the disc is truncated at 1/3 
of the binary separation, and that there is no replenishment of 
the circumstellar discs from circumbinary material. The vertical 
dashed line shows the separation below which the discs suffer no 
mass loss. 



the fraction of those stars which fail to accrete a significant 
new disc following their last close fiyby. 



in the initial conditions for discs than it does about the 
internal physics governing their evolution. 



3.3 Evolution of the accretion rate of CTTS 

With the caveats noted above, the evolution of the disc frac- 
tion in Taurus is consistent with a theoretical model in which 
differences in disc lifetime are ascribed to differences in the 
initial disc mass. A la dispersion of around 0.5 in the log 
of the initial disc mass is required to reproduce typical disc 
lifetimes of 1-10 Myr. Figure |7| shows how this model com- 
pares with observations of the accretion rate of CTTS in 
Taurus. We plot estimates of M by Gullbring et al. (1998) 
(as tabulated in Hartmann et al. 1998) against our own re- 
vised determinations of the stellar ages. The vertical error 
bars show Hartmann et al.'s (1998) best estimate of the er- 
ror in determining the accretion rate, which is 0.6 in log(M). 
Against this data, we plot theoretical curves for the (3 = 3/2 
disc model. We show the predicted evolution of M for the 
central estimate of the initial disc mass of 0.1 Mq, and for 
the models with initial disc masses differing by ±lo\ 

From the Figure, it is clear that the ±lcr curves bracket 
the observed range of CTTS accretion rates at a given age. 
The dispersion in initial disc mass required to reproduce the 
observed range in disc lifetimes also produces roughly the 
correct spread in accretion rates. Since we can reproduce the 
same result with other choices of (3, our interpretation is that 
the Hartmann et al. (1998) plot of M against t, shown in a 
modified form as Figure |7| tells us more about the variation 



4 THE INFLUENCE OF BINARY 
COMPANIONS 

Most T Tauri stars are found in binary systems (Leinert et 
al. 1993; Chez, Neugebauer & Matthews 1993; Simon et al. 
1995), which can affect the evolution of circumstellar discs in 
several ways. First, gravitational torques from the compan- 
ion act to truncate the disc (Paczynski 1977), and prevent 
it from expanding freely. Second, if the binary is sufficiently 
close, it will limit the reservoir of gas available to be ac- 
creted. Both of these effects would be expected to reduce 
the lifetime of discs in binary systems as compared to single 
stars. A substantially reduced disc lifetime in binaries is not, 
however, observed (Simon & Prato 1995; Bouvier, Rigaut & 
Nadeau 1997; Chez, White & Simon 1997; White & Chez 
2001). 

Fig. |H| shows how binaries alter the predicted disc life- 
time in our models. We have assumed that the only effect 
of the binary on the initial surface density distribution is to 
truncate the disc at a radius r± = d/3, where d is the binary 
separation. The initial accretion rate in the binary models 
is therefore the same as for the isolated disc models. At rt, 
we model the torques from the companion using a v r = 
boundary condition. Parameters for the disc wind and the 
disc viscosity remain unchanged. 

By construction, the disc lifetime for sufficiently wide 
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binaries tends to the observed value of « 6 Myr. At smaller 
separations, all three disc models display similar behaviour. 
The disc lifetime initially increases, reaching a peak at a sep- 
aration of 15 au to 100 au, depending upon the model. This 
increase occurs because the disc, which is prevented from ex- 
panding freely, has less area from which mass can be lost in 
a wind. This compensates for, and initially overwhelms, the 
tendency for smaller disc to have reduced viscous timescales. 
As a result, only very close binaries (d less than 10 au to 
40 au, depending upon the model) are predicted to lead to 
a reduction in the disc lifetime as compared to the lifetime 
of an isolated disc. 

In addition to asking whether the disc lifetime around 
the primary is a function of binary separation, we can also 
ask whether differential evolution of the primary and sec- 
ondary disc is expected to yield 'mixed' binary systems that 
pair CTTS and WTTS. The naive expectation is that the 
lower mass component of the binary, whose Roche lobe en- 
closes a smaller disc, should lose its disc and become a 
WTTS first. Observations, however, show that mixed pairs 
are rather rare. Using a sample of close binaries (most of 
which had projected separations between 15 au and 100 au) 
Hartigan & Kenyon (2002) found that only 25 percent (4 
out of 16 systems) of binaries containing at least one CTTS 
were mixed systems. Combining these results with earlier 
studies (Brandner & Zinnecker 1997; Prato & Simon 1997; 
Duchene et al. 1999), which obtained similarly small mixed 
fractions, there is no evidence for a trend in the fraction 
of mixed systems with separation. Comparing the observa- 
tions to disc models that do not include winds, this is only 
consistent with theoretical expectations if the initial condi- 
tions yield more massive discs around the lower mass star 
(Armitage, Clarke & Tout 1999). 

If we assume that the initial disc accretion rate, and 
the mass loss rate, are the same for both stars in a binary, 
then the disc models presented in this paper predict a low 
level of mixed WTTS / CTTS binaries. For a mass ratio 
q — 0.5, the maximum radii of discs surrounding the pri- 
mary and secondary differ by around 50 percent (Papaloizou 
& Pringle 1977). From Figure |S] we then estimate that an 
initially doubly strong (CTTS + CTTS) binary with sepa- 
ration d > 10 au will appear as a mixed binary for at most 
10% (if /3 = 3/2) to 20% (if /? = 1) of the time. A larger 
fraction (up to 50%) of mixed binaries is predicted for the 
(9 = 1/2 disc model. As was the case with the disc lifetime, 
therefore, only very close binaries or binaries with extreme 
mass ratios are expected to show significant differential disc 
evolution. 



5 SUMMARY 

We have shown that protoplanetary disc models which in- 
clude mass loss from the disc at large radius are generally 
consistent with observations of accretion in T Tauri stars. In 
particular, we can simultaneously reproduce both the rapid 
disc dispersal required by observations (Simon & Prato 1995; 
Wolk & Walter 1996; Duvert et al. 2000), and the near- 
independence of the disc lifetime on the binary environment 
(White & Chez 2001). The evolution of the disc fraction 
in young clusters, which we have studied in detail for Tau- 
rus, can also be modelled within this framework, provided 



that we allow for a spread in initial disc masses around a 
central value of ~ 0.1 Mq. We believe that the required 
level of variation in the initial conditions of discs is a plau- 
sible outcome of the turbulent conditions that give rise to 
star formation within clusters. Simulations of clustered star 
formation (e.g. Bate, Bonnell & Bromm 2003) show that in- 
teractions between small groups of forming stars and their 
discs are almost inevitable. In such a dynamic environment, 
we expect large star to star variations in the initial mass 
and angular momentum content of circumstellar material. 
In fact, a significant fraction of the stars formed in the Bate 
et al. (2003) simulation suffer encounters with other stars 
that are close enough to truncate circumstellar discs to sizes 
of au or smaller. Such interactions provide a possible origin 
for the otherwise puzzling observation of a population of ap- 
parently very young Weak lined T Tauri stars seen in the 
Taurus region. 
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